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Abstract 



A Monte Carlo simulation of the 0(4) A</> 4 theory in the broken 
phase is performed on a hypercubic lattice in search of an / = 1, 
J = 1 resonance. The region of the cutoff theory where the interaction 
is strong is investigated since it is there that a resonance would be 
expected to have a better chance to form. In that region the presence 
of an / = 1, J = 1 resonance with mass below the cutoff is excluded. 
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The question of whether or not an / = 1, J = 1 resonance exists in 
the broken phase of the four component A0 4 theory in four dimensions is in 
fact an old one. In the sixties, a quite phenomenological theory describing 
the low energy pion processes was developed by Gell-Mann and Levy [p]]. 
The theory is described by a chiral Lagrangian, and involves the three pion 
fields 7T 1 , 7r 2 , 7r 3 , the scalar field a, the two nucleons, their interactions, and 
the dimensionful pion decay constant f w = 90 MeV. The a, although it 
has never really been seen, may exist in nature as a broad and quite heavy 
(~ 700 MeV) isospin 1=0, spin J=0 resonance. In that setting, it was natural 
to ask whether the presence of the p resonance was a consequence of the pion 
interactions or of some higher energy physics. Since the theory of Gell-Mann 
and Levy becomes an 0(4) A0 4 theory if the nucleon fields are neglected, 
this question reduces to whether or not the 0(4) A0 4 theory can sustain an 
/ = 1, J = 1 resonance in the broken phase. Earlier attempts P] to answer 
this question used analytical approximations such as Pade approximants and 
found that an / = 1, J — 1 resonance is present in the theory. It is not clear, 
however, whether or not the presence of this resonance is an artifact of the 
approximations used and therefore the question of the existence of an / = 1, 
J = 1 resonance in the theory has yet to receive a conclusive answer. 

Today, in a very different setting, this same question is of interest again. 
The scalar sector of the Minimal Standard Model is also a four component 
A(/> 4 theory in the broken symmetry phase. The equivalent of the a resonance 
is the Higgs particle, the three pions are the Goldstone bosons, the pion decay 
constant is the weak scale fa = 246 GeV, and the scattering of longitudinally 
polarized vector bosons behaves exactly the same way as a scaled up version 
of ix-ix scattering (equivalence theorem ||). It is possible that the Higgs, 
like the a, is quite heavy and broad and may avoid detection at SSC. On 
the other hand, if the four component A0 4 theory does indeed contain an 
I = 1, J — 1 resonance, then since the p resonance, as it appears in nature, 
is quite strong, its equivalent in the scalar sector of the Minimal Standard 
Model may have a good chance to be detected at SSC. Therefore, it becomes 
very important to know if another "signature" of the scalar sector, besides 
the Higgs resonance, is awaiting discovery at SSC. 

A direct Monte Carlo simulation was decided to be the best way to shed 
new light on this old question (a leading order large N calculation was not 
able to produce an I = 1, J = 1 resonance). The 0(4) A0 4 theory was 
simulated on the lattice in the A — > oo limit. In that limit, the theory has 
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the strongest interactions and a resonance probably has a better chance to 
form. The model has the lattice action 



s = E \- k E + $ A-a) + a(44 - 1) 2 + (i) 

and was simulated on a hypercubic lattice A of spatial extension L and 
time extension L t using an incomplete heat bath algorithm on the CM-2 
machine at SCRI. 

Since, on a finite lattice, the direction of the symmetry breaking changes 
from configuration to configuration, the same approach as in || was used to 
disentangle the massive scalar field a from the Goldstone modes. For each 
configuration, the global 0(4) coordinate system was rotated so that its first 
axis would be parallel to the direction of the symmetry breaking. Because 
there are many ways to perform this rotation, a "simple" one was consistently 
used throughout the simulation. The field, expressed in the new coordinate 
system, is § x = (a x , 7r*, 7T 2 ., ir%), with a being the massive scalar field and 
7T 1 , 7r 2 , 7r 3 the three massless pion fields. 

To measure the mass m a of the a field, the time slice connected correlation 
function Co,o(t) =< O a (0), O a {t) > c of the zero 3-momentum operator 

Oa(t) ( 2 ) 

where A t is the three dimensional time slice at Euclidean time t, was fitted 
to 



C(t) = a cosh 



+ b (3) 



by a three parameter correlated x 2_ fit with M = m a . For each fit, the errors 
were obtained by varying x 2 by 1. 

To measure the energy of the lowest laying state in the / = 1, J = 1 chan- 
nel, an operator carrying these quantum numbers needs to be constructed. 
The simplest such operator is: 

O c>m (t) = E K^l+m e ab C (4) 

where summation over repeated indices is assumed, a, b, c are the isospin 
indices, e a b c is the totally antisymmetric tensor, rh is the m'th Euclidean 
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unit vector of the time slice At, and m G [1,2,3] is the z-component spin 
index. Unfortunately, the time slice connected correlation function of this 
operator gives a very weak signal. To get a better signal an operator that 
extends over several lattice spacings needs to be constructed using a trial 
wave function for the two-pion state. The "bag" and "bound state" meson 
wave functions were considered |J. The latter gave a better signal with an 
affordable cost in computer time and was therefore used for the simulation. 
Using as a trial wave function for the two pions at relative position R, the 
hydrogen wave function ^ n =2,i=i,m(R), an J = 1, J = 1 operator with total 
3-momentum zero was constructed: 

Oc,m(t) = E I E 1^1 exp(-|i?|/2ao) Y hm (6, ^K^+^abc \ (5) 
sGA ' [rgb J 

where B is a three-dimensional cubic box centered at the origin and contained 
in A t , cio is the "Bohr radius" in lattice units, 9 and <fi are the azimuthal 
and polar angles of R, and Yi jm is the / = 1 spherical harmonic (up to a 
multiplicative normalization constant). The parameter ao can be given any 
value. A very large ao will cause the exponential to decrease very slowly 
and then the sum over R will have to be carried out over a box B as large 
as A t . Since this can increase the computer time significantly, a smaller 
cto has to be used so that the size of the box that contains the important 
contribution from the exponential can be made smaller. However, ao cannot 
be made very small because the signal becomes weaker as ao decreases. An 
optimal choice of ao and B was found to be ao = 2, and B extending from 
—3 to +3 in each of the three directions. The operator Cjm (t) couples to 
the I = 1, J = 1 states. The energy of the lowest laying state in this channel 
can be found by looking at the time slice connected correlation function 
C c , m {t) =< O C)m (0)O C)m (t)* > c of this operator. 

The simulation was done on an L = 8, L t = 16 and L = 16, L t = 16 
lattice and for three values of the hopping parameter n = 0.305, 0.310, 0.330. 
These values were chosen so that a wide range of m a will be covered (they 
also coincide with some of the values used in ||). The expectation values of 
Co t o(t), C\ y i{t) = | J2 c ,m C c ,m(t), and also of the pion zero 3-momentum time 
slice correlations were measured. The pions will not concern us here except 
to mention that they are the massless Goldstone modes and found to behave 
as in [||. 
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K 


L 


sweeps 


measurements 


m a 


0.305 


16 


4xl0 5 


1.6xl0 4 


0.225(4) 


0.310 


16 


1.6xl0 5 


1.6xl0 4 


0.424(6) 


0.330 


16 


0.8xl0 5 


1.6xl0 4 


0.83(1) 


0.305 


8 


17.5xl0 5 


7xl0 4 


0.343(2) 


0.310 


8 


7xl0 5 


7xl0 4 


0.468(3) 


0.330 


8 


1x10 s 


2xl0 4 


0.91(2) 



Table 1: Number of measurements and m a . 



K 


L 


E 


fit-range 


X 2 /d.o.f. 


0.305 


16 


0.94(1) 


2-8 


8.8 


0.305 


16 


0.76(3) 


3-8 


1.8 


0.310 


16 


0.94(1) 


2-8 


0.7 


0.310 


16 


0.91(3) 


3-8 


0.7 


0.330 


16 


0.94(1) 


2-8 


4.2 


0.330 


16 


0.87(2) 


3-8 


2.4 


0.305 


8 


1.549(5) 


1-5 


1.7 


0.310 


8 


1.531(6) 


1-5 


1.7 


0.330 


8 


1.501(9) 


1-4 


0.4 



Table 2: Energy E in the 1=1, J = 1 channel. 

The number of sweeps, the number of measurements of C 0)0 (t) and Ci^t), 
and the the values of m a , for the various lattice sizes and k's, are given in 
table 1. The results for L = 16 are in good agreement with the results of |J. 

Ci i(t) is plotted versus t in figure 1 for the L = 16 lattice and the three 
values of k. Ci t i(t) was fitted with the expression in equation [3] with b = 
for a few different ranges of t, and the resulting energies E = M, together 
with the x 2 P er degree of freedom for each fit are given in table 2 for the 
L = 16 and L = 8 lattices. The dotted line in figure 1 is the fit for k = 0.310, 
L = 16, and 2 < t < 8 and is presented to give the reader a feeling of the 
quality of the fits. 
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Figure 1: < C\^{t) > for the L = 16, L t = 16 lattice. 

The effective energy E e s(t), obtained by fitting Cx,i(t) to the expression 
in equation^, with 6 = for the two time slices at t-1 and t, is plotted versus 
t and for the three values of k in figures 2 (L — 16) and 3 (L = 8). The 
values of t omitted from those plots had an E e g(t) with very large error. 

In a two-pion I = 1, J = 1 state with zero total 3-momentum, the lowest 
3-momentum a pion can have has one component equal to ^ and two equal 
to 0. The next one has two components equal to ^ and one equal to 0. The 
energy spectrum in the I — 1, J — 1 channel is therefore expected to contain 
levels with energies close to the energies of these states, but slightly different 
because of the interaction. For the L — 16 lattice, these levels have energies 
Eq ~ 0.78 and E± ~ 1.09 respectively, and are denoted by the dotted lines 



6 



2.0 



1.5 



1.0 



0.5 



o k = 0.305 
□ k = 0.310 
o k = 0.330 




0.0 



2 



4 



6 



Figure 2: E cS (t) for the L = 16, L t = 16 lattice. 

in figure 2. From this figure, it is clear that the observed energy levels are 
very close to the levels of two free pions. In fact, for smaller t, the levels are 
close to E\, and for larger t they are close to Eq. Because the free two-pion 
levels are not very well separated at L = 16, the observed levels are probably 
a mixture of the two lowest ones. In that sense, the energies in table 2 for 
the L = 16 lattice are probably a mixture as well. The fact that the observed 
levels correspond to a two-pion state and not to a resonance is also greatly 
supported by the fact that these levels change only slightly from k = 0.305 to 
k = 0.330. After all, in that range m a varies from 0.225 to 0.91. Therefore, 
if a resonance is present it must have energy larger than 0.78 and is either 
too heavy (for example, heavier than 1.09) to be observed with this method, 
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Figure 3: E cS (t) for the L = 8, L t = 16 lattice. 

or is "hiding" between 0.78 and 1.09. That the latter is not the case can be 
seen by looking at the energy levels obtained for the L = 8 lattice. From 
figure 3 it is clear that there are no energy levels below ps 1.4. In fact, since 
the two lowest free two-pion states are well separated in this lattice, the 
lowest energy level is clearly visible in figure 3. It is true that the limited 
statistics give E eS (t) only up to t = 4, but because the lowest level is well 
separated from the next one, the correlated x 2 fit gives a good estimate of 
the energy of this level. The numbers given in table 2 are indeed very close 
to the lowest free two-pion energy E ~ 1.50 and do not seem to change 
much for the different values of k. From this analysis it is evident that if an 
/ = 1, J = 1 resonance exists for k > 0.305, it is unphysical since it must 
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have energy in lattice units greater than 1 (physical energy above the cutoff). 

For the theory of Gell-Mann and Levy of low energy pion processes the 
k > 0.305 region corresponds to where the a particle mass is greater than 
approximately 180MeV (or equivalently the cutoff is less than approximately 
1.3GeV).[] Since m a is expected to be much larger than 180MeV, it is clear 
from the results that the existence of the p resonance in nature cannot be 
accounted for by this low energy theory alone (in contrast to 0). 

For the scalar sector of the Minimal Standard Model the k > 0.305 region 
corresponds to where the Higgs mass is greater than approximately 500Ge\^ 
(or equivalently the cutoff is less than approximately 3.5TeV). The Higgs 
mass is of course not known but its upper bound is placed at around Q50GeV 
(hypercubic lattice triviality bound). Thus the existence of an I = 1, J = 1 
resonance can be excluded with confidence for values of the Higgs mass above 
« 500GeV. For « < 0.305 (Higgs mass < 500GeV) the strength of the 
interaction becomes weaker and hence it is safe to say that if a resonance 
could not form for k > 0.305, where the interaction is stronger, it is unlikely 
that it will in this region either. For this reason it was not deemed necessary 
to investigate the k < 0.305 region. A numerical simulation for k < 0.305 not 
only is not necessary, but it would also be very costly since larger lattices will 
have to be used (the correlation length becomes larger than approximately 



It must be emphasized that these conclusions are valid only within the 
realm of the scalar sector of the Minimal Standard Model. It is of course 
still possible that the physics that enters at higher energies may be able to 
produce such a resonance in very much the same way the physical p particle 
owes its existence to QCD. This resonance if it exists due to some higher 
energy theory it would have energy determined by that theory. In fact, it is 
possible that the energy of this resonance is determining the cutoff energy of 
the Minimal Standard Model. 
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1 The connection with the physical scale f v (or fa) is made through the renormalized 
coupling constant gn = 3m^/f%. The value of gu for this value of k was taken from 0. 
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